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From a consideration of high temperature series expansions in ferromagnets and in spin glasses, 
we propose an extended scaling scaling scheme involving a set of scaling formulae which express to 
leading order the temperature (T) and the system size (L) dependences of thermodynamic observ- 
ables over a much wider range of T than the corresponding one in the conventional scaling scheme. 
The extended scaling, illustrated by data on the canonical 2d ferromagnet and on the 3d bimodal 
Ising spin glass, leads to consistency for the estimates of critical parameters obtained from scaling 
analyses for different observables. 
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Critical divergences of thcrmodynamical quantities 
F(T) at continuous phase transitions are convention- 
ally quoted in terms of the normalized scaling variable 
t={T- Tc)/Tc, i.e. in the form 



F{T)^Af[{T~T,)/T, 



(1) 



where Tc, p and Ap are the transition temperature, 
the critical exponent, and the critical amplitude, respec- 
tively. There exist associated finite size scaling (FSS) 
rules. It has been clearly underlined (see e.g. [l|) that 
this representation is valid only in the immediate vicin- 
ity of Tc, which is a very restrictive condition both for 
numerical simulations and experiments. In particular for 
the analysis of finite size numerical data on complex sys- 
tems such as Ising spin glasses (ISGs) , where simulations 
have many intrinsic limitations, it is difficult to analyze 
data while complying strictly to this condition. 

Other scaling variables can be used. In all modern 
theoretical and numerical analyses on ferromagnets, e.g., 
Ref. I^ldi &s well as in some experimental analyses, e.g., 
Ref. dill, the scaling variable t = 1- (3/ = {T -Tc)/T 
is used instead of t. Although no general rationalization 
seems to have been published yet explaining why one 
scaling variable should be chosen rather than another, in 
this Letter, we propose a coherent scaling scheme for crit- 
ically divergent quantities derived from a systematic con- 
sideration of high temperature series expansions (HTSE) 
which naturally leads us to use the variable r. Beside 
this, the HTSE analysis leads us to properly define sin- 
gular terms of interest in such a way that they themselves 
reproduce appropriate temperature dependence at high- 
est temperatures, i.e., in the limit /3 — > 0. We call this 
the extended scaling scheme which we demonstrate below 
to be quite powerful at temperatures close to Tc where in 
practice one makes critical analyses for divergent quan- 
tities. 



For a ferromagnet and in particular an Ising ferromag- 
net (IF), our extended scaling scheme is explicitly de- 
scribed as follows. 

(i) We use t as the normalized scaling variable, and write 
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for the reduced susceptibility x following the standard 
definition without a prefactor /?, which is in fact consis- 
tent with the idea of our extended scaling, 
(ii) Defining the second moment correlation length ^ 
through M2_= J2r''"'^ i^oSr) — 2fixC^j with d the spatial 
dimension|3, we write down ^ as 



m 



(3) 



(iii) Using this form of ^(/3), we rewrite the FSS ansatz, 
F{L,(i)^LP/-F[Lim],''^s 

F{L, P) ^ [L/(3^'^r/^:F [(L//?i/2)i/^(i _ ^/^^)] . (4) 

For the 2d-l¥, where the confluent corrections to scal- 
ing are known to be zero 0, S 13 ; we demonstrate that 
the above critical expressions with the known exact val- 
ues of the critical parameters reproduce x(/3) and ^(/3) to 
a good approximation right up to high temperatures, and 
that our FSS form Eq. holds to a high approximation 
over a very much wider range of L and T than the conven- 
tional one. For the ISC, 0^ replaces /? throughout in (i), 
(ii) and (iii), giving an entirely novel set of expressions. 
The use of our FSS form appropriately modified for the 
ISG resolves a longstanding puzzle in the ISO critical 
analysis, i.e., published estimates for the critical expo- 
nent V through X scaling and through ^ (or the Binder 
parameter) scaling differ by a large factor j^, . We em- 
phasize here that Eqs. ©^101) above are leading order 
expressions of the divergent quantities. This does not 
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mean that we neglect even the confluent corrections to 
scahng, but that fits of the numerical data examined be- 
low to our extended scaling scheme only with the leading 
expressions are quite satisfactory. We will discuss sepa- 
rately ^3 that analyses of published high precision data 
on canonical ferromagnets using the present extended 
scaling scheme actually give estimates of the confluent 
correction terms which improve considerably over those 
from standard analyses. 

In standard spin 1/2 ferromagnets the HTSE for the 
susceptibility xiP) is written as 



xiP) = 1 + + 02/3' + a3/3^ + • • • 



(5) 



where /3 = J IknT with the coupling constant J and fee 
set to unity |3l. lllUl2l |. The asymptotic form of its factors 
a„ is eventually dominated by the closest singularity to 
the origin (Darboux's First Theorem which in the 
simplest case is the physical singularity, i.e.. 



Leading term 
From Ref. [2] 
FromRef. [12] 




t={T-Tc)/Tc 



Eq.(2) 
From Ref. [2] 
FromRef. [12] 
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FIG. 1: Susceptibility x of the 2d-IF as a function of log(t) 
(a) and log(r) (b). 
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with being the inverse critical temperature. One of 
the techniques to relate these two expressions is the ra- 
tio method, in which the recurrence relation a„/a„-i = 
(l//3c)(l + (7 - ^)/n) for large n is used Q- It is there- 
fore natural to adopt t ~ 1 — (3 //3c as the scaling variable 
in critical analyses based on the HTSE theory. 

The HTSE for the second moment /i2(/3) introduced 
above is of the form Q 



(7) 



It diverges at as [T — T^]"'''"*"^'^-'. Then, invoking again 
Darboux's theorem to link the series within the brackets 
[■ • ■] to the critical divergence, the appropriate extended 
scaling form can be written as 



-(7+21^) 



(8) 



Combined this with Eq. © for xiP), Eq. ® for ^(/3) 
is derived. With this expression for ^ the ESS form be- 
comes Eq. Q. At the limit P ^ 0, Eqs. Q, ® and 
© have the same /? dependence as the leading terms of 
the corresponding HTSE. This implies that they merge 
smoothly to the analytic corrections to scaling to yield 
the proper expressions at highest temperatures. Similar 
expressions are expected also for confluent corrections to 
scaling if they exist. 

We exhibit in Figs, ^a) and^^b) log-log plots of the 
susceptibility x of the canonical 2d-lF in the thermo- 
dynamic hmit, plotted against \og{t) and log(r), respec- 
tively. The data points in the flgure are the high-precision 
results of the critical and HTSE analyses. The 
line in Figs, ^a) and (b) is the power-law expression of 
X as a function of t and r, respectively, with 7 = 7/4 
and the critical amplitude = 0. 962581. ..0- By the 



T scaling, all the data points nearly up to r = 1, i.e., to 
almost T = 00, lie on the scaling expression of Eq. (0) 
within the accuracy of the figure. This result reflects the 
weakness of corrections to scaling in this system as men- 
tioned above. On the other hand, the deviation of the 
true X data from the t scaling line, x — AyJ:^'^ , is sig- 
nificant already at, say, t ~ 0.2 (or T ~ 1.2Tc). Thus 
if the scaling variable t rather than r were used for a 
critical analysis on this system, there would appear to be 
very strong "correction" terms, which would pollute the 
evaluation of the 7 value. 

Figure 13 shows the ESS plots for the 2d-\F susceptibil- 
ity, the standard one in the inset and the extended one 
in the main frame. It is clear that the standard form, as 




FIG. 2: An extended FSS plot for the 2d-IF susceptibility 
X- The points are data obtained by our Monte Carlo (MC) 
simulation at temperatures 0.46 < T/Tc < 1.86. The inset 
presents a standard FSS plot for the same x- For both plots 
the exact scaling parameters are used. 
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^{T,L)/L 

FIG. 3: An extended (conventional) FSS plot for the 2d-IF 
susceptibility x normalized by {L\/tY~'^ {L^~'^) as a function 
ofC(r,L)/L. 



expected, gives acceptable scaling only extremely close to 
Tc, while the extended form gives high quaUty scaling for 
all temperatures above Tc examined. In Fig. |31we show 
the conventional and extended FSS plots of x{T,L) as a 
function of £,{T,L)/L. Note that, since the value Tc is 
not involved in this analysis at all (also in Fig. Elbelow), 
one can judge straightforwardly our proposal (iii) by this 
comparison. The consequence is that the extended FSS 
plot is definitely better than the conventional one. We 
thus conclude that, at least for the 2d-IF, our extended 
scaling scheme does indeed work much better than the 
conventional one. 

Let us next discuss an extension of our extended scal- 
ing scheme to the Edwards- Anderson 3d-ISG model hav- 



ing a symmetric interaction distribution with zero mean. 
Because of the symmetry, as stated by Daboul et al 0| , 
only even powers of /3 enter into the HTSE for thermo- 
dynamic quantities, and the HTSE for the reduced spin 
glass (SG) susceptibility XSG (the ordinary one multi- 
plied T^) is of the form 



Xsg(/3) = 1 + ci(/32) + c2(/32)2 + csipy + 



(9) 



Hence once again invoking the Darboux theorem, but this 
time with (/3//3c)^ replacing /3//3c, or r' = 1 — W/ ffr)"^. 
we adopt the following scaling form for xsoiP) 



Xsg(/3) ^ A^sG (1 - iP/Pc)') 



2\ -1 



(10) 



The //2 in spin glasses, again due to the symmetry, can be 
expressed as even powers of /3 starting from the term, 
though the coefficients have not been explicitly evaluated 
yet. Hence the scaling form for ^(/3) in spin glasses can 
be taken to be 



(11) 



Then the extended FSS for ISGs can be written as 
F{L, (3) ^ [L/PY/"T [{L/fif'^ (1 - ifi/Pcf)] . (12) 

In contrast to the 2d-l¥, analytical theories are very 
limited for ISGs. Daboul et al [l5j were able to make ac- 
curate estimates of (3c and 7 of the models by the HTSE 
method but only in dimension 4 and above. We therefore 
compare our extended r scaling scheme with the conven- 
tional one with t variable, without introducing correc- 
tions to scaling in either analysis. The numerical data 
used are obtained on the 3d-ISG system with bimodal 
interactions by the exchange MC method. 
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FIG. 4: An extended FSS plot of ^{L,T)/L of 3d-ISG. The 
straight line with slope 2.72 {— v) represents an expected 
asymptotic form of the scaling function. The inset shows the 
FSS plot of the Binder parameter. 
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FIG. 5: An extended FSS plot of xsG- The inset is the cor- 
responding conventional FSS plot. 
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FIG. 6: An extended (conventional) FSS plot for the 3d-ISG 
susceptibility XSG normalized by (LT)^"'' (L^"'') as a func- 
tion of c(r,L)/L. 

In Fig.Elwc show an extended FSS plot for the corre- 
lation length £^{L,T)/L as a function of the scaling vari- 
able X = It'HTL)^^" based on Eq. In this plot 
we fix Tc = 1.11, which is the optimal value from our 
analyses and is consistent with that of Ref. [3, and ad- 
just v to obtain the best scaling fit. We end up with 
f = 2.72(8). The fit is surprisingly good for all the data 
with L's indicated in the figure and at T from O.SlTc to 
7.2Tc. The slope of the straight line in the range |a;| 3> 1 
is 2.72 (= v). Note that data points on the line arc not 
only those at sufficiently high temperatures with £, ~ f3 
from Eq. ^11^ but also those in the critical range with 

Figure [S] shows an extended FSS plot for XSG- This 
plot is obtained by fixing Tc = 1.11 and v = 2.72 and by 
adjusting rj to give rj = —0.40(4). We obtain quite sat- 
isfactory scaling for all our data. In contrast, a conven- 
tional FSS plot using the variable t as shown in the inset 
is rather poor except for the immediate vicinity of x = 0. 
The fit yields a small apparent value of = 1.47(3)simi- 
larly to the previous conventional FSS analyses |3, 1^ . In 
Fig. El we demonstrate conventional and extended scal- 
ing plots for xsG versus ^. The comparison of the two 
implies that our scaling scheme with (ii) and (iii) is defi- 
nitely more appropriate also for the 3(i-ISG. 

One can remark that our extended scaling scheme with 
only the leading term for each divergent observable gives 
high quality fit over the entire range of L and T exam- 
ined. Another important result here for the SG study 
is that the extended FSS analyses on ^, g and xsG with 
the implicit assumption that corrections to scaling are 
weak yield a unique critical parameter set. This is in 
sharp contrast to standard FSS methods for which the 
estimate for v obtained from x scaling with the same as- 
sumption is considerably smaller than that from S. or g 
scaling d, 0- In this context, we note that numerical 
data on the same 3rf-ISG model have also been analysed 



using t scalin g to gether with strong correction to scal- 
ing terms flH Il7l|. We have checked that our data can 
equally well be analysed using a very similar method to 
that of Ref. and the excellent fit to the data from 
the extended scaling without corrections to scaling is ob- 
tained with two fewer fitting parameters. This does not 
imply that corrections to scaling are absent, but just that 
their influence on the fits examined here is rather weak. 

In conclusion, by considering the intrinsic structural 
form of high temperature series developments, we have 
proposed an extended scaling scheme with appropriate 
scaling expressions for thermodynamic observablcs in fer- 
romagnets and in spin glasses, and have demonstrated 
the results which support it strongly. One of them is the 
direct comparison of our extended scaling on x vs. ^ with 
that of the conventional one. Another is the result that, 
within our scheme, the leading order critical power-law 
expressions with a coherent set of critical parameters re- 
main good approximations to the true behavior over a 
much wider temperature range than with the standard t 
scaling. From these results we consider that our extended 
scaling scheme with the variable r is more fundamental 
than the conventional t scaling. 
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